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I. INTRODUCTION 

Braneworld models of the universe — in which the observable universe is a four dimen- 
sional timelike hypersurface (brane) embedded in a higher dimensional (bulk) space-time — 
have attracted much recent attention. This is partly due to the fact that Superstring/M- 
theory seems to require the existence of extra dimensions and the braneworld approach 
may be one way of reconciling our 3+1 dimensional universe with these higher dimensional 
theories l|, [j, [Sj . 

Another reason for the current popularity of the braneworld construct is due to the fact 
that brane cosmology is usually accompanied by new features and is therefore, in principle, 
falsifiable 0, Q]. The simplest Randall-Sundrum (RS) braneworld, for instance, gives rise 
to an evolutionary equation for the brane which differs from standard general relativity at 
early times This leads to several interesting consequences. For instance, the very early 
universe expands as oc p, instead of the more familar H oc in standard cosmology. 
The changed expansion rate causes a scalar field to experience greater damping, which, in 
turn, allows Inflation to occur for a broader class of initial conditions and potentials j^. If 



ly bounces and avoids the initial 
(The behaviour of anisotropies 



the fifth dimension is timelike then the universe generica 
big bang singularity which plagues standard cosmology 

in the RS scenario can also be very different from that in standard general relativity 

A complementary approach to braneworld cosmology pioneered by the DGP model |9(], 
allows the universe to accelerate at late times thus providing a geometrical answer to the 
riddle posed by dark e nerg y. Models which unify the RS and DGP approaches also lead to 
several new features j^, [lo{ |. For instance (i) the phenomenon of dark energy can be transient 
so that the universe accelerates for a while before settling back into matter dominated 
expansion, (ii) the effective equation of state of dark energy can be phantom-like {wcs < 



— 1), (iii) new cosmological singularities can arise in such models Such alternative 

cosmological models provide reasonable fits to the current cosmological data [l^ . 

In this paper we address the issue of cosmological evolution on a brane in a theory of 
gravity whose action includes, in addition to the familiar Einstein term, a Gauss-Bonnet 



contribution. Gauss-Bonnet terms arise naturally in superstring theories 13|] and their 
cosmological effects have been discussed in several papers lisl. [igI. The present paper 
deals with this issue in greater generality, we examine both cases: when the bulk dimension 
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is spacelike as well as timelike. We also a develop a new pictorial method of analysis which 
provides qualitative insights into the evolution of the universe in this potentially important 
new model of gravity. 

II. BASIC EQUATIONS 



We begin with the following n-dimensional (n > 5) action: 

1 



S 



-9 



{R-2A + aLGB) 



(1) 



where R is the n-dimensional Ricci scalar, A is the n- dimensional cosmological constant, and 
Kn '■= y/SirGn, where G„ is the n-dimensional gravitational constant. The Gauss-Bonnet 
term Lgb is a combination of the Ricci scalar, the Ricci tensor and the Riemann tensor 



^GB 



R' 



(2) 



The constant a in ([T]) is the coupling constant of the Gauss-Bonnet term and for a — ^ our 



model reduces to the familiar Randall-Sundrum model [3|. The action ([T]) can be obtained in 
the low-energy limit of heterotic superstring theory 13|], in which case a can be regarded as 



the inverse string tension and is positive-definite. We, therefore, assume a > throughout 
this paper. (We shall explicitely be assuming n > 5 since for n < 4 the Gauss-Bonnet term 
is a topological invariant and does not contribute to the field equations.) 
The gravitational equations which result from the action ([1]) are 



+ aH^"^ + A^'^^ = 



where 
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(3) 

(4) 

(5) 



A. Bulk solution 



The n-dimensional vacuum solution can be obtained as a product manifold ^ x 
^i^]-^ ^]^Q }i]2e element 



ds^ = —hir)dt^ + e——— + r'^'Jijdx^dx^ , 
h{r) 



(6) 
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where is an (n — 2)-dimensional space of constant curvature with unit metric 7jj. In 

the equations which follow, k denotes the curvature of and takes the values 1 (positive 
curvature), (zero curvature), and —1 (negative curvature). The value of the constant e 
determines whether the (bulk) fifth dimension is spacelike (e = 1) or timelike (e = — 1). In 
the former case, is a Lorenzian manifold, whereas in the latter case, it is a Euclidean 
manifold. 

The basic equations of the theory under consideration are 







2a{n - 3){n - A)h - e\r^ + 2ak{n - 3)(n - 4) 



dr'^ 



+2(n-3)r 



2a{n - 4)(n - 5)/i - elr'^ + 2ak{n - 4)(?i - 5) 



dh 
dr 



+2ar^{n - 3){n - 4) 



fdhV- 



+a{n - 3){n - 4)(r2 - 5){n - 6)h^ 



-e(n-3)(n-4)/i 



+ 2ak{n-5){n-6) 



-2Ar^ + k{n - 3)(n - 4)r2 + ak^{n - 3)(n - 4)(n - 5)(n - 6) 



(7) 



2a{n - 3)(n - A)h - elr'^ + 2ak{n - 3)(n - 4) 



{n — 2)r 



+a{n - 2)(n - 3)(n - 4)(n - 5)h^ 



dh 
dr 



-e{n - 2){n - 3)h 



+ 2ak{n - 4)(n - 5) 



-2Ar^ + k{n - 2){n - 3)r^ + ak^{n - 2){n - 3){n - 4)(n - 5) 



(8) 



where the former is the component of Eq. ([3]), while the latter is the (t,t) or (r, r) 
component acting as a constraint. The general solution of these equations is 

^2 



h{r) = ek + 



ST 



afi 8{n-3){n-A) 



\ r"-i {n-l){n-2) 



aA 



(9) 



2{n - 3){n - A)a 

where /x is a constant. Our solution for h{r) has two branches, which correspond to the two 
signs in front of the square root in Eq. iQ. We call the family with the minus (plus) sign 
the minus-branch (plus-branch) solution. 

• For e = 1, the minus-branch solution has the general relativistic limit as a 0, 
while there is no general relativistic limit for the plus-branch solution. (The global 



structures of this solution were presented in 



For e = —1, the plus-branch solution has the general relativistic limit as a 0, while 
the minus-branch solution does not. 
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Hereafter, we shall be considering a five- dimensional bulk spacetime, for which the metric 
(EI) reduces to 



dr 
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dsl = g^^dx^'dx" = -h{r)dr + e—— + [rfx + fk{x) {dO' + sin^ ed(j)')\ , (10) 

ii[r) 



Mr) = .fc + ^(.T^l + ^ + ^«A), (11) 

where /o(x) = X, hix) = sinx, /-i(x) = sinhx, and e = ±1. 

In this spacetime, there are two classes of singularities when fi ^ 0. One is the central 
singularity at r = and the other is the branch singularity at r = rb := [—afi/{l + 
4aA/3)]^''^ > 0, when the term inside the square-root in Eq. (fTTj) vanishes. The branch 
singularity exists if fi is negative, or if 1 + 4a A/3 < for positive /i. 



B. Priedmann equation on the brane 

The position of the three-brane is described by the functions r = a{T) and t = T{t) 
parametrized by the proper time r on the brane. The tangent vector to the brane is written 
as 

d • d , d 
^ dx>^ dt ~^ ' dr ' ^ ^ 

where a dot denotes the differentiation with respect to r. The normalization condition 
u,,u^ = — 1 leads to 

-2 

1 = h{a)f' - , (13) 

and the induced metric of the three-brane cjab is given by 

dsl = Qabdy^dy^ = -dr'^ + a(r)^ [dx^ + fkixfidO^ + sin^ ed(j)'^)] . (14) 

The unit normal 1-form to the three-brane is given by 

n^dx^ = ddt — Tdr , (15) 

where n^n^ = and n^n^ = 1/^ satisfied. 

The extrinsic curvature of the three-brane is obtained from Kab '■= n^.^e'^e'^, where : = 
dx^/dy"'. We have 

e1dy'' = idT, eldy" = ddr , e^dy" = d'.dy^ , (16) 
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and 



Then, we obtain the non-zero component of K'^i, as 

K\ = — ^ I a + — 1 , ICi = 
hT\ 2eJ ' ' 

where a prime denotes differentiation with respect to a. 



hT ■ 




ea 



The junction condition at the brane is given by 
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{n) 



(18) 



(19) 



where 



Jab := \ [2KKadK\ + K^fR'^f Kab~2KadK''^ Kfb- K^Kab 



P, 



adbf 



Radbf + 2ha[fRb]d + 2hd[bRf]a + Rha[bhf]d - 



(20) 
(21) 



The energy-momentum tensor T°'b on the brane is given by 



T^-b = diag{- p,p,p,p) + diag(-a, -cr, -a, -a, ) 



(22) 



where p and p are the energy density and pressure of a perfect fluid on the three-brane, and 
the constant a is the brane tension. We have introduced the notation 



[X]±:=X+-X- 



(23) 



where is the quantity X evaluated either on the + or — side of the brane, and Padbf is 
the divergence-free part of the Riemann tensor, i.e.. 



DaP"'dbj — . 



(24) 



where Da is the covariant derivative on the brane. 

From the (r, r) and (i,i) components of Eq. (fT9l) and Eq. (fT3|) . we obtain 



K 



^ ^ 4a / 3A:-g/i(a) ^ 



(25) 



where H := a/ a. Here, we have assumed Z2-symmetry of reflection with respect to the 
brane. This generalized Friedmann equation reduces to that obtained by Davis [l^ for 
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Differentiating Eq. (125|) with respect to r and using Eq. (fT3|) and the (r, r) and 
components of Eq. f|T9|) . we obtain 

p = -3i/(p + p), (26) 

which is the energy-conservation equation on the three-brane. Let us assume that the perfect 
ffuid on the three-brane obeys 

p=(7-l)p, (27) 

where we assume that the equation of state of matter on the brane hes within the Zeldovich 
interval < 7 < 2 due to the dominant energy condition (equivalently, —l<w<l, where 
w : = p/ p = 1 + '-f is the equation of state). From Eq. ( 1261) . we then obtain 

where po is a positive constant, so that p is a monotonically decreasing function of a for 
7 > 0. 

C. The Randall Sundrum limit 

In this paper, we shall consider only those solutions of fllip and 0251) which possess the 
general-relativistic limit since other solutions may describe physically inadmissible evolution 
of our brane. The minus- and plus-branch solutions in (fTTl) have the general relativistic limits 
for e = 1 and —1 in (!25|) . respectively. 

As mentioned earlier, the action ([T]) contains the Randall-Sundrum model as a subclass. 
Setting a ^ in Eq. ( l25l) . one gets the generalised Randall-Sundrum (RS) model 

'iiP + ^y = ^(k-^- W) + sH' . (29) 



36 a2 V 8a2 6 

From this equation and from Eq. fl28l) . we obtain 



^f= = J:f^ + -y-^H-^4A- (30) 



3Qe Va37 J Sa^ 6 

The Randall-Sundrum model corresponds to e = 1, while the dual model with e = — 1 was 
discussed in Permitted values of the expansion factor must clearly satisfy > 0. An 
interesting consequence of (l30l) is the possibility of singularity-free solutions when e = —1 







8 



III. PICTORIAL ANALYSIS OF COSMOLOGICAL EVOLUTION 

We saw in the previous section that the evolution equation for the Gauss-Bonnet brane 
can be quite comphcated and, therefore, difficult to analyze analytically. In this section, 
we present a general method of analysis which allows one to study pictorially the behaviour 
arising from the generic cosmological equation (l25l) . 

We notice that equation fl25l) . describing the cosmological evolution of the Gauss-Bonnet 
brane, always has the form of a cubic curve in the {H^, p^ot) plcine: 

Cpl, = {A±H^){B + Hf , (31) 

where ptot '■= p + cr , A and B are functions of a, C is a positive constant and the ± sign 
corresponds to e = ±1. The value of cosmological constant a can be positive, negative or 
zero. The right-hand side of equation ( pTl) has exactly three real zeros in H^, two of which 
coincide, namely, {H^)i = l^A, and {H'^)^^ = —B. Only part of this cubic curve lies in the 
physical domain > 0, p^^^ > 0, and it is in this domain that the evolution of the brane 
takes place. Consequently, the evolution of our brane-universe can be pictured as a point 
moving along this cubic curve in the physical domain p^^^ > 0, > 0. 

This pictorial representation is very useful in appreciating the full gamut of possibilities 
for cosmic evolution of this brane. For comparison, it is helpful to note that cosmological 
evolution in general relativity (GR) is described by 



k 



= Ptot - ^ , (32) 



where we have set the proportionality term SttG/S to unity. Equation (l32l) describes a 
quadratic curve in the (if^, p1^^) plane. Another example is the Randall-Sundrum brane, 
which, for every value of a, is described by a straight line in the (i?^, Pto^) plane: 



^' = — + + (33) 



where e = ±1, and we have set the term ^5/86 in fl30l) to unity. 

Before commencing our discussion on the subtleties of cosmological evolution on the 
Gauss-Bonnet brane, it will be helpful to first consider the different evolutionary possibilities 
in a spatially flat universe described by the more familiar general-relativistic equation (l32ll 
with k = 0, where a acts as a cosmological constant. In this case the expansion of the 



9 



universe can proceed in three distinct ways, corresponding to the cases cr > 0, o" < and 
(7 = 0. All three possibilities correspond to motion along the quadratic curve in Fig. [H 

Notice that expansion along the entire curve from the top (early times) to the origin (late 
times) takes place only if a < 0. In the absence of a cosmological constant (a = 0) the origin 
(if^, p^oJ = (0, 0) is reached at r — > cxd. In contrast, for a < 0, the origin is reached in a 
finite interval of time when the matter density has dropped to p = At this point H = 0, 
in other words expansion ceases and the universe begins to contract. Evolution thereafter 
proceeds upward along the same curve — in reverse fashion. 

Finally, if a > 0, evolution does not proceed all the way to the origin but terminates 
at some point D along the curve. At this point, p ^ so that ptot = cr and = a. 
The end point D of evolution corresponds to the universe's asymptotic approach towards 
de Sitter space. (This, for instance, would be the case for a spatially flat ACDM universe 
which accelerates at late times.) 



+ 




h2 



FIG. 1: The evolution of a spatially flat FRW universe in GR proceeds along this curve. The 
downward arrow indicates expansion while the upward arrow indicates contraction. The latter is 
only possible if o" < 0. For o" > the expansion of the universe terminates at the point D at which 
p = 0. At this point, the universe expands exponentially. For cr = 0, the origin {H = 0, p = 0) 
marks the end point of evolution. The scale of the x and y-axis is arbitrary. 
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A. Spacelike extra dimension (e = 1) 

Let us now discuss the evolution on the Gauss-Bonnet brane in greater detail. For a 
spacelike extra dimension, the cosmological equation ( l25l) has the form 

Cpl, = {A + H^){B + Hf , (34) 

where 
and 

in general, are functions of the scale factor a. 

As mentioned earlier, equation (l34l) has the form of a cubic curve in the {H^, p^ot) plane. 
The two signs in fl36|) correspond to the two different ways of embedding the brane in the 
bulk space. In this paper we only consider the upper sign, which has the GR limit. 

As discussed in the previous section, the evolution of the braneworld is described by a 
point moving along the cubic curve in the {H"^, Ptot) plane, in the physical domain > 0, 
Ptot — 0) with the parameters of the cubic curve simultaneously changing with time due to 
the dependence of A and B on the scale factor (see below). The evolution can proceed in 
three distinct ways which are summarized below. All three cases correspond to 5 > in 
equation fl36l) . and the first two also have A < 0. 

1. The behaviour of the universe is shown in the left panel of Fig. [2l The point P 
corresponds to = —A. During the course of expansion, the motion along the curve 
is initially downwards from the initial Big Bang (BB) singularity towards P. However, 
for P to be reachable in a finite time interval the brane tension a must be negative 
since only then is (p + a)^ = permitted. The point P marks a turning point for the 
evolution along the curve: after this point, the energy density of the universe keeps 
decreasing while the quantity (p + cr)^ is increasing. In the case p = 0, we also have 
H = at the point P. In this case, the Hubble parameter passes through an inflection 
point at P. Since 

- = H + H\ (37) 
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it follows that a = aH^ > when p = \a\. In other words, > for some length 
of time during the upward motion along the curve away from P. Thus the universe 
accelerates at late times. Note that H = —AnGp in ACDM and if > is usually 
associated with a phantom equation of state w < — 1 in dark-energy models. (This 
qualitative behaviour will take place also for sufficiently small values of afi/a^ reached 
in the course of expansion in the neighbourhood of the point P, which will make A 
almost constant in time.) 

The growth in H, however, cannot continue indefinitely since p ^ eventually, and 
(p + (t)^ (corresponding to the point D), which implies — > const. This means 

that the universe approaches a de Sitter-like attractor (D) at very late times. We 
therefore conclude that our Gauss-Bonnet brane can display phantom-like features 
and super-accelerate at late times, before approaching w — — 1 in the distant future. 
Note that the big-rip future singularity (at which H oo) is absent in this case, 
which is one of the appealing features of this scenario. (Other braneworld models with 



this property have been discussed in |lOl ]. 



2. For cr > 0, the point P can never be reached, and expansion proceeds along BB D, 
culminating in de Sitter-like asymptotic expansion at D. 



(a) (b) 

FIG. 2: Spacelike extra dimension: A < 0, B > (left) and A > 0, B > (right) in The 
point P is the turning point, and the point E is the point of recollapse. 
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3. The case with B > and A > is shown in the right panel of Fig. [2l The point E is 
the point of recoUapse. At this point the (spatially fiat) universe ceases to expand and 
begins to contract. The point E is reachable either if the brane tension is negative, or 
if it is positive with the value of o"^ lying below the point E. 

One should note that, theoretically, the scale-factor dependent parameter A{a) can change 
sign during the course of evolution, so that the curve along which the evolution takes place 
can continuously evolve from that in the right panel of Fig. [2] to that in its left panel, and 
vice versa. This introduces an obvious modification to the description of the evolution, 
which does not change in any significant way. 

The complete set of figures showing the {H^, p^J plane are shown in Fig. [6] of the 
Appendix. 

B. Timelike extra dimension (e = —1) 

Also in this case, there is only one branch of the generic cosmological equation ( l25l) having 
the GR limit which has the form (!3T]) . namely, 

Cpl, = {A-H'){B + H^f , (38) 

where 
and 

A:=^(jl + ^ + ^aA-l] , 5 := ^ f 1 + l./^T^^A^U - + - . (40) 
4a \^ V 3 J ' 4a \^ 2 V 3 j 8a 2 ^ ' 

Clearly, the theory makes sense only for A > (for which the branch singularity does not 
appear), hence, also i? > 0. A typical graph illustrating the case < i? < 2 A (equivalently 
l/4a < A) is shown in Fig. [3](a). The point P corresponds to = A. The graph 
corresponding to i? > 2A, which is equivalent to 

^<i' " ^IH-f + 5"A<2, (41) 

is shown in Fig. [3](b). 
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The end points E and P in Fig. [3]^a) are the reverse points of evolution. The point S is 



the position of a sudden 'quiescent' singularity of the type described in [11|, Il9|. Indeed, the 
evolution of the universe cannot be continued beyond this point because the quantity p^^^ 
should change in the same direction (decrease), which is physically impossible. Note that 
the value of H is finite and nonzero at this point, while H is divergent. (This can easily be 
seen by writing d{H^)/dT = d{H'^)/dptot ■ dptot/d^ where c/ptot/T = —<iHp and noting that 
d{H^)/dp tot — ^ oo at S*.) The Kretschmann invariant on the brane K := Rabcd.R"'^'^'^ is given 
by 

K = 12[H^ + {H + H'^f] , (42) 

and diverges as the quiescent singularity is approached. 

Consider now in more detail the evolution of the Gauss-Bonnet brane suggested by 
Fig. [3]^a). The Big Bang singularity which featured prominently in Fig. [2] has effectively 
been replaced by the sudden singularity S. The following four possibilities for evolution 
immediately suggest themselves: 

1. Expansion commences at S and proceeds to which marks a turning point at which 
H = 0. Thereafter, the universe ceases to expand and begins to contract. The con- 
tracting trajectory ends at S. The sudden singularity as S marks both the beginning 
and end of evolution in this scenario. (The possibility that quantum effects might 
modify cosmolo gica l evolution in the vicinity of such a singularity has been discussed 



in l20l; see also 



2l|.) 



2. The universe contracts from the singularity at S until it reaches E, where it bounces, 
then expands back to S. In this case, the brane tension must be negative {a < 0) since 
that is a necessary condition for moving along the trajectory SE during contraction. 
In the vicinity of p = \cr\, the map (p + cr)^ — p is bivalued (see Fig. H]), which allows 
(p + cr)^ to increase both when p increases as well as decreases. This ambiguity is 
responsible for the two possibilities discussed above. 

3. The trajectory S ^ P describes a super-accelerating universe expanding from the 
singularity S, since it suggests that increases while p decreases. (In fact, H oo 
at the point S.) If p = 0, then H > throughout this phase, and it is unlikely 
that SP in this case can describe the real universe. If p > 0, then, in the course of 
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the evolution, super-acceleration may be replaced by the "usual" acceleration. If the 
brane tension is negative, then the point P is reached, after which the evolution turns 
back to the PS path. Then, depending on the value of the brane tension, it either 
reaches the singularity S again or asymptotically approaches the de Sitter state at an 
intermediate point between S and P. 

4. For a brane with negative tension, one also has the time-reversal of the previous 
case, which describes a universe contracting either from the de Sitter state at an 
intermediate point between S and P or from the singularity S*, proceeding to P and 
then to singularity 5". 




(a) (b) 
FIG. 3: Gauss-Bonnet brane with timelike extra dimension: 

(a) < i3 < 2A. The points E (bounce or recollapse) and P are the turning points of the evolution, 
while the point S corresponds to a sudden or 'quiescent' singularity. 

(b) B > 2A. The bouncing scenario (which requires o" < 0) describes a brane contracting from a 
de Sitter- like initial stage at P to Q at which p = |(7| and ptot := p + o" = 0. Further contraction 
takes the universe from Q to E, and along this segment both p as well as p + a increase. At E, the 
density of the universe has reached its maximum value while the Hubble parameter has declined 
to zero. The universe therefore bounces at E, then re-expands and evolves in reverse fashion along 
E ^ Q ^ P. Note that P marks the beginning and end point of evolution. 

We remember that during the evolution the parameters A and B describing the cubic 
curve change their values, and it may happen that the curve changes its shape during the 
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u < 




P 



FIG. 4: For negative values of the brane tension (cr < 0) two values of the matter density p map 
onto a single value of ptot = /? + cr, as illustrated in this figure. 

evolution, that some critical points leave the physical domain H"^ > 0, p^^^ > or, on the 
contrary, enter this domain. All such possibilities are quite easy to investigate case by case, 
but we will not do this in this paper. 

The complete set of figures showing the [H"^, Ptot) plane are shown in Fig. [7] of the 
Appendix. 

C. Bouncing Braneworld 

In order to address the issue of a bounce in the Gauss-Bonnet brane in more detail, let us 
first consider this issue within the context of the Randall-Sundrum model (which presents 
a limiting case of our braneworld). As mentioned earlier, cosmological evolution of the RS 
brane is described by equation fl33p which represents a straight line in the {H^, pl^^ plane. 
We show this line in the left panel of Fig. Elfor a time-like extra dimension [e = —1). 

One can see a close qualitative similarity between the curve on the left panel of this figure 
and the curve in Fig. [3](b). The bouncing scenario in figures [3](b) and [5] proceeds as follows: 
the universe begins to contract from the point P at which p = and H = const. In other 
words, both the starting point and end point of evolution correspond to de Sitter space. 
An increase in the value of the matter density brings us to the point Q at which p = \a\ 
and ptot := p + cr = 0. (Note that cr < is a prerequisite of this model, since otherwise 
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(a) (b) 

FIG. 5: The RS brane with a timehke extra dimension. Note that the left panel shows ptot = 
as a function of H'^ while the right panel shows as a function of . Since the map (p + ci)^ — > p 
is bivalued (see Fig. S]), it is easier to discern the salient features of the bounce in the right panel 
than in the left ! 

the point Q cannot be reached.) The universe contracts further from Q to E and in this 
segment both p as well as p + cr increase. At E the density has reached a finite maximum 
value while the Hubble parameter has declined to zero. The universe therefore bounces at 
E, then re-expands and evolves in reverse fashion along E ^ Q ^ P. 

IV. DISCUSSION 

Braneworld cosmology has attracted considerable interest during the past decade. This 
is partly due to the fact that such models may play an important role in the low energy 
limit of M-theory/string theory. Another reason for the growing interest in brane dynamics 
is associated with the new features which some of these models possess and which, in turn, 
can lead to new cosmological predictions and scenario's. Our attempt in this paper has 
been to develop a completely general qualitative approach to determine the salient features 
of a brane embedded in a five dimensional bulk and evolving according to the precepts 
of Einstein-Gauss-Bonnet gravity. For this purpose we show that the 3+1 dimensional 
equations of motion of several popular cosmological models can be depicted as simple curves 
in the {H"^, {p + o")^) plane. (Here H is the Hubble parameter, p the density and a the brane 
tension.) 
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For instance, the spatially flat FRW universe in GR has the form of a quadratic curve 
while the Randall-Sundrum model describes a straight line in the {H"^, (p + cr)^) plane. The 
Gauss-Bonnet brane, on the other hand, describes a cubic curve in the {H"^, {p + cr)^) plane 
— see equations ( l33l) . ( l32i) and ( 1311) . This pictorial depiction of dynamics permits us to 
discover the salient features of cosmic evolution very simply. Applying this approach to the 
Gauss-Bonnet brane we discover the following interesting properties: 

1. For a flnite region in parameter space the Gauss-Bonnet brane accelerates at late 
times. Acceleration can be phantom-like {w < —1) but does not lead to the eventual 
destruction of the universe in a big-rip future singularity. Instead, at very late times 
the expansion of the universe approaches de Sitter space and becomes exponential 
(i.e. w —1). (The possibility that the current expansion of the universe may 



be phantom-like has evoked much recent interest and discussion; see 



10 



22| for a 



non-exhaustive list of papers discussing this issue and 23j for a summary of recent 
observational results.) 

2. The expansion of the universe may commence from or terminate in a 'sudden' quiescent 
singularity, at which the Hubble parameter and the density of matter remain flnite 
but H diverges. 

3. The universe can evade the initial big bang singularity and bounce. (This possibility 
is realized if the flfth dimension is timelike.) 

Whether any of these properties of the Gauss-Bonnet braneworld is realised in practice 
is currently an open question which can be answered by: (i) a deeper understanding of 
the embedding of this cosmology within a more fundamental theoretical framework, (ii) a 
comparison with observations. 
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V. APPENDIX 

Figures [6] and [7] describe the evolution of our brane with a spacehke and a timehke extra 
dimension, respectively. These figures supplement those appearing in the main body of the 
paper. 

(a) xP^, (b) 




FIG. 6: The (i?^, Ptot) plane is shown for the GB brane with a spacelike extra dimension for the 
following values of the parameters A, B in ()36p with the upper sign: (a) A < —B/2, (b) A = —B/2, 
(c) -B/2 < A < 0, (d) A = 0, (e) < A ^ B , and (i) A = B (at the branch singularity). Here, 
xi := —{B + 2A)/3 and X2 ■= —B. Note that the region with < or p'^^^ < is nonphysical. 
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(c) 




FIG. 7: The (H^, 

Ptot) plane is shown for the GB brane with a timeUke extra dimension for the 
following values of the parameters A, B in (00]): (a) B > 2 A, (b) B = 2 A, and (c) < i? < 2 A. 
The strange singularity characterized by yl = and B = 3/(8a) appears in the extremal case of 
(a). Here, x-^ := {2A — B)/3 and X4 := —B. Note that the region with < or p^^^ < is 
nonphysical. 



